Suppressing nano-scale stick-slip motion by feedback 
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When a micro cantilever with a nano-scale tip is manipulated on a substrate with atomic-scale roughness, the 
periodic lateral frictional force and stochastic fluctuations may induce stick-slip motion of the cantilever tip, 
which greatly decreases the precision of the nano manipulation. This unwanted motion cannot be reduced by 
open-loop control especially when there exist parameter uncertainties in the system model, and thus needs 
to introduce feedback control. However, real-time feedback cannot be realized by the existing virtual reality 
virtual feedback techniques based on the position sensing capacity of the atomic force microscopy (AFM). To 
solve this problem, we propose a new method to design real-time feedback control based on the force sensing 
approach to compensate for the disturbances and thus reduce the stick-slip motion of the cantilever tip. 
Theoretical analysis and numerical simulations show that the controlled motion of the cantilever tip tracks 
the desired trajectory with much higher precision. Further investigation shows that our proposal is robust 
under various parameter uncertainties. Our study opens up new perspectives of real-time nano manipulation. 



I. INTRODUCTION 

One of the central problems in nano science and tech- 
nology^ is the realization of high-precision nano manip- 
ulation, e.g., pushing, pulling, rotating, rolling, and cut- 
ting nano-scale objects. In the widely applied AFM ex- 
periments 2 13 , a large amount of frontier progresses have 
been achieved about nano manipulation. However, the 
theoretical analysis for such system is mainly focused 
on the static force analysis within a Newton-mechanical 
framework 4 ^"— that describes the macroscopic system. 
Recent experiments show that the force analysis in the 
nano scale is not the same as that at the macroscopic 
scale. Typically, nano-scale frictio n 14 ! 15 is logarithmi- 
cally dependent on the velocity^ - —, which is quite dif- 
ferent from the macroscopic sliding friction (the friction is 
independent on the velocity) and the traditional viscous- 
type friction (the friction is linearly proportional to the 
velocity). Additionally, the atomic scale periodic struc- 
ture on the substrate has to be seriously considered in the 
friction analysis that is usually done with the continuous 
mechanics in macroscopic syste m 19 ! 20 . 

The periodic lateral force induced by the atomic peri- 
odic structure on the surface of the substrate may lead to 
stick-slip motions 21 23 of the cantilever tip of the AFM, 
which deteriorates the precision of the nano manipula- 
tion. In the literatures^ 5 the cantilever of the AFM is 
usually modelled as a soft spring, and the effective nano 
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friction force is described by a sinusoidal lateral force. 
Temperature-dependent white noises are introduced to 
represent the stochastic fluctuations in the lateral force 
induced by the thermal motion of the substrate atoms. 
Such a dynamical model predicts a motion transition of 
the cantilever tip of the AFM from the continuous slid- 
ing to a stick-slip mode by varying the ratio between the 
amplitude of the sinusoidal lateral force and the stiffness 
of the cantilever tip. 

To improve the precision in nano manipulation, many 
strategies 24 30 have been proposed to control the motions 
of the nano objects under nano frictions. For example, in 
the literature^ 4 , the authors designed a non-Lipschitzion 
control function under which one can push the nano sam- 
ple to asymptotically track the target velocity. This pro- 
posal is efficient and robust, but the natural fluctuation 
cannot be removed. To overcome this difficulty, more 
complex control function was designed based on the Lya- 
punov theory in Refs. 29,30 . However, the control designs 
in these proposals require the knowledge of the exact po- 
sition of the sample during the nano manipulation. Such 
schemes are uneasy to be realized with the present ex- 
perimental techniques, e.g., the haptic sensing and the 
virtual reality visual feedback techniques ^ 10 ! 31 , for the 
inability of simultaneous position sensing and manipula- 
tion processing by the AFM. In fact, in these techniques, 
one has to stop the nano manipulation process and scan 
the surface of the substrate by the cantilever tip of the 
AFM to position the sample on the substrate, after which 
the next step of nano manipulation can go on. To solve 
this problem, in this paper, we propose a feedback con- 
trol strategy based on the real-time signal sampled from 
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the force sensor of the AFM, which can be conditionally 
done without pending the nano manipulation. The sig- 
nal is used to estimate the position of the cantilever tip 
of the AFM, with which we can design feedback control 
to reduce the stick-slip motion of the cantilever tip and 
thus improve the manipulation precision. 

The paper is organized as follows. Section HI1 describes 
the dynamical model, following which the stick-slip mo- 
tions of the cantilever tip under open-loop control are 
presented in Sec. Illli Section HVl is devoted to the design 
of the real-time feedback control to reduce the stick-slip 
motion of the cantilever tip. Section [V] discusses the ro- 
bustness of our method against various parameter uncer- 
tainties. Conclusions and forecast of the future work are 
given in Sec. I VI I 



II. MODELLING OF THE ONE-DIMENSIONAL NANO 
MANIPULATION SYSTEM 

We first present the model used to describe the one- 
dimensional nano manipulations such as pushing a nano 
sample or etching the surface of the substrate to draw 
desired pattern. In this model, the cantilever of the AFM 
is taken as a spring with an effective stiffness k c . Thus, 
according to the Hooke's law, the lateral force imposed 
by the cantilever tip is: 



k c (u - x), 



(1) 



where x and u are the relative positions of the cantilever 
tip with and without deformation on the platform of the 
nano manipulation, respectively (see Fig. [1]). 

In our method, u is the control parameter to be de- 
signed. The mechanism of control is shown in Fig. QJa). 
As shown in Fig.QJa), in the designed nano manipulation 
system, the cantilever of the AFM is fixed, while the plat- 
form of the nano manipulation moves. In such a system, 
the motion of the platform of the nano manipulation is 
controlled by a piezoelectric transducer (PZT). By ad- 
justing the voltage Vpzt added on the PZT via a control 
circuit, the deformation of the PZT can be controlled, by 
which the relative position u between the cantilever of the 
AFM and the platform is tunable. Typically, the func- 
tional relationship between the deformation of the PZT 
and the voltage Vpzt shows hysteresis, creep, and struc- 
tural vibration behaviors and thus is nonlinear. However, 
such a nonlinear characteristic response of PZT can be 
compensated by auxiliary control devices. Our recent 
theoretical and experimental study 32 shows that a lin- 
ear functional relationship between Vpzt and u can be 
obtained by introducing optimal design of feedforward 
controller by the Prandtl-Ishlinskii model 33 . Based on 
this study, in order to simplify our discussions, we do 
not go into details about the relationship between Vpzt 
and u, but simply take the deformation of the PZT repre- 
sented by u as the control parameter. Additionally, in the 
moving frame of the platform, the control system shown 



in Fig. QJa) is equivalent to that shown in Fig. QJb), in 
which the platform is fixed and the cantilever moves. 



Origin of the ^a*^ 
coordinate ""*—*"' " : 




FIG. 1. (color online) Schematic diagrams of the controlled nano ma- 
nipulation system: (a) the control system in experiments in which the 
cantilever of the AFM is fixed and the platform moves; (b) the equiva- 
lent control system in which the platform is fixed while the cantilever 
moves. The two control systems (a) and (b) are equivalent if we fix the 
origin of the coordinate with the moving platform. "PZT" denotes the 
piezoelectric transducer. The deformation of the PZT is controlled by 
the voltage Vpzt added on the PZT. Thus, the relative position u can 
be controlled as what we want by a control circuit. 



When the cantilever tip moves on the substrate, the 
nano friction may occur on the contact area. In recent 
experiments 16 , the average frictional force is observed to 
be logarithmically dependent on the velocity of the can- 
tilever tip v = dx/dt, i.e., 



■/o 



F h In 



(2) 



where F/ , Ff t , v\ are constant parameters. There also 
exists periodic lateral force induced by the surface poten- 
tial of the substrate, and stochastic noises induced by the 
thermal motions of the substrate atoms. For simplicity, 
we only consider the fundamental-frequency component 
of the periodic lateral force, and omit the correlation ef- 
fects of the thermal motions of the substrate atoms. Ad- 
ditionally, to simplify our discussions, the origin of co- 
ordinate is chosen such that the periodic lateral force is 
zero at the origin. Thus, the resulting modification of 
the lateral force can be expressed as: 



SF f = F h sin (27^) - £(t), 



(3) 



where Ff 2 is the amplitude of the periodic lateral force; a 
is the lattice constant of the substrate; and £(t) is a white 
noise such that E(£(t)) = 0, E(£(t)€(t')) = D^6(t - t') : 
with Dg. being the strength. E (•) denotes the ensemble 
average over the white noise. 

From Eqs. $1]), ((2]), and ((3]), the equations of the system 
can be expressed in the Ito notation as: 



dx 
mdv 



vdt, 
k c (u 



x)dt - 
F /2 sin(2^) 



Ffo 
dt- 
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D 6 dW, 



(4) 
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where m is the effective mass of the cantilever tip; and 
dW = W (t + dt) — W (t) is the increment of the Wiener 
process W (t) = J* £ (r) dr satisfying: 

E(dW) = 0, (dW) 2 = dt. (5) 

III. STICK-SLIP MOTION UNDER OPEN-LOOP 
CONTROL 

In order to implement nano manipulation by the can- 
tilever tip of the AFM with a constant velocity v* , i.e., 
to control the position of the cantilever tip such that 
x = v*t, a simple strategy is to set u = v*t. However, 
this would be ineffective due to the existence of the pe- 
riodic lateral force and the stochastic fluctuations espe- 
cially when the roughness of the substrate is relatively 
large compared with the stiffness of the cantilever of the 
AFM. In such case, existing studie o 19 i 2Q have predicted 
stick-slip motions of the cantilever tip, which greatly re- 
duce the precision of the nano manipulation. The lateral 
pushing force imposed by the cantilever tip oscillates un- 
der the stick-slip motion, which can be so large that the 
fragile sample and substrate are damaged, or the can- 
tilever tip of the AFM slide over them. 

The drawbacks of the open-loop constant control can 
be seen from the following numerical examples. The sys- 
tem parameters are chosen as^: 

F fo = 10 nN, F h = 1 nN, a = 0.25 nm, 
vi = l nm/s, F h = either 0.25 nN or 15 nN, 
= 3 nm/s, yfDl 0.1 nN • s" 1 / 2 , 
k c = 1 N/m, m = 5 x 10 _n kg, (6) 

under which the simulation results of the position x, ve- 
locity and the lateral force F p imposed by the can- 
tilever tip are shown in Fig. [2] 

As shown in Fig. [2j when the amplitude of the peri- 
odic lateral force Ff 2 is relatively small such that Ff 2 /a 
is comparable with k c (the case with Ff 2 = 0.25 nN), 
which is the case when the cantilever tip moves on an in- 
commensurate substrate^, there exists no stick-slip mo- 
tion and the cantilever tip moves with a constant velocity 
after a transient process (see the green triangle curves). 
However, when the amplitude of the periodic lateral force 
is relatively large such that Ff 2 /a ^> k c (the case with 
Ff 2 = 15 nN), which may be valid when the cantilever 
tip moves on a commensurate substrate, stick-slip mo- 
tions can be observed (see the blue solid curves). Such 
stick-slip motions can be explained by the elastic insta- 
bility predicted by the Frenkel-Kontorova modeU^. In 
this case, the cantilever tip moves only when the pushing 
force F p exceeds the critical value, i.e., the peak value of 
the sawtoothlike trajectory shown in Fig. [2fc). It is also 
shown in Fig. [2fa) that no matter whether the stick-slip 
motion occurs or not, there exist tracking errors between 
the controlled trajectories (the green triangle curve and 
the blue solid curve) and the ideal trajectory (the red 




FIG. 2. (color online) Plots of (a) the positions, (b) the velocities of 
the cantilever tip, and (c) the pushing forces imposed by the cantilever 
tip. The red curves with plus signs denote the ideal trajectories for 
which the cantilever tip moves with constant velocity v = v* — 3 
nm/s; the green triangle (blue solid) curves represent the controlled 
trajectories under open- loop control x = v*t, with the amplitudes of 
the periodic lateral force Ff 2 as 0.25 nN (15 nN). 



curve with plus signs), which need to be compensated to 
improve the precision of the nano manipulation. 



IV. REDUCTION OF THE STICK-SLIP MOTION BY 
REAL-TIME FEEDBACK CONTROL 

In this section, we introduce feedback control to re- 
duce this unexpected stick-slip motion under open-loop 
control. The object of our method is to monitor the de- 
formation of the cantilever and thus make the cantilever 
tip move with a given constant velocity. 

To acquire feedback signals, there are two available 
sensing methods supported by the AFM-based nano ma- 
nipulation systems: position sensing and force sensing. 
Position sensing is one of the basic functions of the AFM, 
by which one can easily obtain the tomography of a sur- 
face with nano-scale roughness. However, real-time posi- 
tion sensing is unavailable during the nano manipulation. 
Another sensing capacity of the AFM is the force sensing, 
which can be done by detecting the deformation of the 
cantilever of the AFM. In contrast with position sensing, 
force sensing can be done during the nano manipulation, 
which makes the real-time feedback control possible. The 
main idea of our feedback control method is to estimate 
the position of the cantilever tip by force sensing to adjust 
the control parameter u, and further control the motion of 
the tip. The control process thus can be divided into two 
steps, i.e., a filtering and estimation step and a feedback 
control step, which will be specified below. 
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A. Position estimation by force sensing 

By the force sensing capacity of the AFM, the pushing 
force measured can be expressed as: 



F™ = k c (u-x) + V (t), 



(7) 



where r] (t) is a white noise induced by the measurement 
apparatus, r] (t) satisfies 

E(rj(t)) = 0, E{r 1 {t)r 1 {t')) = D v 5(t-t'), (8) 

where D v is the strength of the noise r\ (t) . 

To reduce the measurement-induced disturbance of 
r] (t) in feedback control, we filter the measured signal 
F™ by a low-pass filter over a time window [t — T,i\. The 
output signal can be expressed as (see, e.g., Refs i 34 i 35 ): 

K (*) = \ f e-V^F? (r) dr, (9) 

where jf t is the damping rate of the low-pass filter. 
Then, the position of the tip x can be estimated from 
F? by: 



1 



(10) 



x(t) = u(t)--F™(t). 

Under the filtering condition 

1ft « 1/T, (11) 



it can be verified that 

x(t)-x (t) = e-v** (x - x ) , 



(12) 



where xq and xq are the initial states of x (t) and x (t) 
respectively (see the derivation in Appendix |A|) . Thus, 
the estimated position x (t) exponentially converges to 
the position of the tip x (t) with the convergence rate 
j ft, and thus x (t) can be taken as a good estimation of 
x (t) when t ^> 1/7/t- 



B. Feedback control design based on the estimated 
position 

Based on the estimated position x given in Eq. (JTJ5J), 
we can design the position-based feedback control u (x) 
to reduce the stick-slip motion of the cantilever tip, which 
is given as follows: 



u (x) = — 



h c x* + F fo + F fl In ( V — j + F h sin ( 2tt^ 



/' 



-k x (x — x*) — ki / (x (r) — x*) dr 



(13) 



where = v*t is the desired motion of the tip; and 
k x , kj > are the control parameters to be determined. 
The term 

k c x* + F fo + F fl ln (v*/v{) + F h sin (2i\xja) 



in the control given by Eq. (fT3|) is introduced to compen- 
sate the mean and the periodic friction forces. The pro- 
portional and integral feedback control terms k x (x — x*) 
and ki J* (x — x*) dr are introduced to speed up the con- 
vergence of the system dynamics to the stationary motion 
and reduce the static error. 




FIG. 3. (color online) Schematic diagram of the feedback control loop, 
where r* = k c x* + F fo + F fl In (v* /vi ) and e = x - x* ; "PZT" , "PEC" , 
"LPF", and "SWG" are the piezoelectric transducer, the photonelectric 
convertor, the low-pass filter, and the sine-wave generator respectively. 
The feedback control circuit is divided into two parts, i.e., a filtering 
circuit and a control circuit The designed feedback control u in Eq. (fl3l) 
can be used to find the voltage Vpzt = Vpzt (u) added on the PZT. 



The schematic diagram of the feedback control pro- 
posed is given in Fig. [3l The deformation of the can- 
tilever tip of the AFM is detected by an optical refract- 
ing system, and then converted into electric signals by 
a photonelectric convertor. The output signal masked 
by the measurement noise r] (t) is fed into a low-pass fil- 
ter followed by a control circuit to generate the feedback 
control signal. The feedback electric signal is used to con- 
trol the motion of the platform by adjusting the voltage 
Vpzt added on the piezoelectric transducer connected to 
the platform. 

Under the filtering condition ([TT]h the weak noise con- 
dition 

<C 2F h (k c + k x ) a 2 /v\ 2F h mv\ (14) 

and choosing the control parameters fcj and k x such that 

fcj < F fl (k c + k x )/mv*, (15) 

the controlled trajectory x (t) tracks the desired trajec- 
tory x* = v*t in average (see Appendix [B]) , i.e., 



lim (E(x(t))-x*) = 0, 



(16) 



where E (•) denotes the ensemble average of the stochas- 
tic signal. This result indicates that the stick-slip motion 
of the cantilever tip can be efficiently suppressed by the 
designed feedback control. 
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To evaluate the magnitude of the stochastic fluctua- 
tion, we further estimate the variances of the position x 
and velocity v of the cantilever tip, which are defined by: 

V x =e((x-E (x)) 2 ) , V V =E (0 - E(v)) 2 ^j . 

With additional calculations, the stationary values of the 
variances V x and V v can be approximately estimated as 
(see the analysis in Appendix IB]) : 



V c 



2F fl (k c + k x y 



v c 



(17) 



The effectiveness of our proposal can be demonstrated 
via numerical examples. Given the system parameters 

F fo = 10 nN, F h = 1 nN, a = 0.25 nm, 
vi = 1 nm/s, Ff 2 = 15 nN, 
= 3 nm/s, 0.1 nN • s 1 / 2 , 

k c = 1 N/m, m = 5 x 10" 11 kg, fe 7 = 1 N/ (m • s), 
fc x = 5N/m, T = 0.1s, 7/t = ls- 1 , (18) 

we compare the motions of the cantilever tip driven by 
the feedback control given in Eq. ([T3]) and by the open- 
loop control u = v*t. Each curve is obtained by averaging 
over 20 sample (stochastic) trajectories. As shown in 
Fig. 2J the stick-slip motion observed under open-loop 
control (the green triangle curves) is greatly reduced by 
the proposed feedback control (the blue solid curves). 
To compare the stochastic fluctuations with the mean 
trajectories, we calculate the square roots of the variances 

1/2 1/2 

Ox — V x and a v = V v . It can be obtained over 20 
stochastic trajectories that = 0.06 nm and = 
0.26 nm/s, which are quite close to the estimated values 
= 0.05 nm, <t£° 0.38 nm/s given by Eq. (H7J) and 
are negligible compared with the average motion (E (x) 
is about tens of nm, and E (v) is about 3 nm/s). 



V. ROBUSTNESS AGAINST PARAMETER 
UNCERTAINTIES 

The system parameters involved in Eq. ([T8]) can be 
identified offline by, e.g., pre-designed nanofriction exper- 
iments. In practical experiments, we also need to con- 
sider the uncertainties in the system parameters which 
may deteriorate the performances of the nano manipu- 
lation. For example, the plastic deformation of the tip 
and the adhesion force between the tip and the substrate 
can lead to small deviation Ak c of the stiffness k c (typ- 
ically Ak c /k c < 10% in the literature^). To reduce the 
effects of uncertainties, we add an integral control term 
—ki Jq (x — x*) dr in Eq. (fT3|) to reduce the static error 
induced by the uncertainties. 

Denote the additive uncertainties in the system param- 
eters k c , F fo , F fl , Ff 2 , v u and a by Afc c , AF /o , AF h , 
AFf 2 , A^i, Aa, and assume that there exists a phase 
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FIG. 4. (color online) Plots of the motions of the cantilever tip with 
(a) the position x, (b) the velocity v, and (c) the lateral force F p im- 
posed by the cantilever tip; (d) shows the trajectories in presence of 
parameter uncertainties given in Eq. (|21f) . The red curves with plus 
signs are the ideal trajectories with v = v* — 3 nm/s; the green curves 
with triangle signs represent the trajectories driven by the open-loop 
control u = v*t; and the blue solid curves denote the mean trajectories 
under the feedback control given in Eq. J 131) . The mean trajectories in 
the plot are obtained over 20 stochastic trajectories. 



offset A(/> in the sine function in Eq. (J3j). With the anal- 
ysis given in Appendix [Cj the static tracking error can 
be controlled to zero, i.e., limt^oo (E (x) — x*) = 0, if the 
uncertainties of the system parameters are not too large 
to satisfy 

A, = \Ak c \ + 2^ | Aa| + 2zl \ A F h \ <k c + k x , 

\AF fl \<F fl , (19) 

and the gain of the integrator ki is chosen such that 

kj < (F fl - \AF fl |) (k c + k x - A x ) /mv*. (20) 

It can be seen that the proportional and integral feed- 
back terms k x (x — x*), ki J* (x — x*) dr in the control 
(fl3j) both contribute to the robustness of our method: (i) 
the integral term is used to reduce the static error; and 
(ii) the proportional term is used to increase the robust- 
ness of our method about the parameter uncertainty (the 
parameter regime given in Eqs. ([19]) and ([20]) is enlarged 
when we increase the control parameter k x ). Given the 
parameter uncertainties: 

AF fo = 1 nN, AF fl = 0.1 nN, Aa = 0.01 nm, 

AF h = 0.1 nN, Ak c = 0.1 N/m, (21) 

Figure Hfd) shows that our method is still valid under the 
given parameter uncertainties (the controlled trajectory, 
i.e., the blue solid curve, matches very well with the ideal 
trajectory, i.e., the red curve with plus signs). 



VI. CONCLUSION 

In summary, we propose a feedback control strategy 
to reduce the stick-slip motion of the cantilever tip in 
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a AFM-based nano manipulation system. The feedback 
control is designed based on the position estimation of 
the cantilever tip obtained by the force sensing capac- 
ity of the AFM. Compared with open-loop control, our 
proposal can greatly reduce the stick-slip motion of the 
cantilever tip. Our method is robust against small un- 
certainties in the system parameters, e.g., the stiffness 
of the cantilever of the AFM, the lattice constant of the 
substrate, and the phase offset in the surface potential of 
the substrate. 

Future study will be focused on extending the method 
to more practical cases. For example, as shown in our 
discussions, our designed feedback control is valid only 
under small uncertainties. More robust design should be 
developed for large uncertainties. Since the uncertainties 
in the model of the system can be reduced to a DC plus 
periodic disturbance, internal model principle is a good 
choice of the control design. Additionally, our control 
design can be naturally extended to the two-dimensional 
nano manipulation systems. More interesting results, 
such as the suppression of the S-shaped motion of the 
cantilever tip, are hopeful to be observed for the two- 
dimensional case. 
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Appendix A: Derivation of Eq. (|12[) 



From Eqs. ([7j) and ([TO]) , we have 



u-x=(F™- V (t)) /k c , 
u-x = F™/k c . 



(Al) 
(A2) 



It can be calculated from Eqs. © and ()A1|) that 
1 r* 

/ e -7/ t (t-r) ( u ( T ) _ a; ( T )) rf r 
T J t-T 

-It* (t-r) v ( r ) dT ^ 




TJt 



/ e -^ f r,(t-f)df 



rj(t — f ) df 



1 J t-T 



dr 



(A3) 



Here, we have used the condition given in Eq. (fTT]h i.e., 
T <C l/7 /t , to obtain e~ 7 ^ w 1 when t G [0,T] and the 
ergodic property-^ of the white noise r] (t) to replace the 
time average of r\ (t) by its ensemble average, i.e., 

1 ( l 

- j V ( V )dr^E( V (t))=0. 

Let us set 

h(t)=u (t) -x(t), 
h(t)=u (t) -x(t), 
then from Eq. (|A3|h i.e., 



h(t) 



t-T 



-r/i(t-r) ft ( r ) dr, 



we have 



jh (t) = ±(h (t) - e~^ T h (t - T)) - Vt h (t) . (A4) 
We assume that g (T) = e lftT h (t + T), then we have 
I (h {t) - e -~ft t T h {t - T)) = 1 (g (0) - g (-T)) , 
which can be replaced by 

g'(T)\ T=0 = jh{t) + 1St h{t) 

under the condition T <C l/jft- Thus, Eq. (|A4|) can be 
reexpressed as: 



From the definition of h and ft, we have 
d 



dt 



(x - x) = - Jf t (x - x) 



(A5) 



(A6) 



which means that x — x — >• when t ^> 1/7/t, i.e., the 
estimated position of the cantilever tip x tracks the actual 
position x in the long time limit. 
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Appendix B: Derivation of Eqs. ([16]) and (fT7) 

By substituting the feedback control law given in 
Eq. (d3j) into Eq. dJ), we have 



dx 
rndv 



vdt, 

[ (x *^ ) (*^ *^ ) 



+F /2 sin 2tt 



in K) 



| In— - In— 



kiq 

d5x = —■jftSxdt, dq = (x — x*) dt, 



dt 



/D ( dW, 



(Bl) 



where Sx = x — x and q = J Q (x — x*) dr. 

By taking the ensemble average and denoting Sx = 
E (x) - £*, Sv = E (v) - u*, Sx = E (Sx), and £ = E (g), 
we can expand the equation to the second-order quadra- 
tures to obtain: 



Sx 
mSv 



Sv, 



2tt ( x'" 
(k c + k x ) Sx + Ff 2 — cos I 27r — ) Sx 



a 



— k x Sx - 

+o 



ELsv + F h ^V x 



2v* 2 

Yl_si_ 

a 4 a 



O 



Sx Sx 
a a 
V 2 

v V 

V* 4 



8tt 3 

3~ 



cos I 2tt — 1 Sx 

a 



O 



Sv^ 

v* 2 



O 



Sx ■ 



-jf t Sx, q = Sx + Sx, 



(B2) 



where V x = E (x — E (x)) 2 and = E (v — (v)) 2 are 
the variances of x and i>; and O (•) denotes the higher- 
order terms. Since the nonlinear equation (jBl|) is lin- 
earized near the origin (0,0,0,0)' in Eq. JB2]), we have 
introduced the Gaussian assumption to omit higher-order 
quadratures of x and v. We can further omit the terms re- 
lated to V x /a 2 and V v /v* 2 in Eq. (|B2|) . In fact, as shown 
in Eq. ([T7J), we have V£° <C a 2 , l^ 00 <C v* 2 under the 
assumption (|T3|) . Thus, we have 14 <C a 2 and K <C v* 2 
for sufficiently long time. Since we are just interested in 
the stationary behaviors of the system dynamics, we can 
omit the terms related to V x /a 2 and V v /v* 2 . 

With the above analysis, we can obtain the lin- 
earization equation of Eq. (jB2j) in the neighborhood of 

(5x,6v,6£,q) = (0,0,0,0)' . It can be easily verified 
that the characteristic equation of the coefficient matrix 
of the linearization equation is 



+ 1ft) 



(k c 



kj_ 

rn 



0. 



Since kj < Ff 1 (k c + k x ) /mv* from Eq. (fT5]h the real 
parts of the eigenvalues of the above equation are all neg- 
ative. It means that the linearization equation is expo- 
nentially stable, and thus the original nonlinear equation 



(|B2|) is asymptotically stable at the origin, which leads 
to the fact that lim^oo (E (x) — x*) = 0. 

To approximately estimate the stationary variances, 
we replace x by x : linearize Eq. (jBl|) , and omit higher- 
order correlation terms. It can be verified that V x , V v , 
and C xv , i.e., the covariance between x and satisfy the 
following equation 







mv* 
1 



2 







H 


IS.) 











(B3) 



To consider the stationary variances, the omission of the 
higher-order correlation terms to obtain Eq. (|B3[) is rea- 
sonable, because x — x — > when t —> oo, and higher- 
order correlation terms are small compared with V x , V v , 
and C xv under the weak noise assumption (fl4|) . From 
Eq. (|B3|) . we can calculate the stationary variances of x 
and v as: 



V c 



D f v* 



2F h (k c + k x y 



v c 



D^v* 



Appendix C: Robustness analysis of our method 

Let us replace the system parameters fc c , Ff , Ff ± , Ff 2 , 
v u and a in Eq. (|B2j by fc c +Afc c , F /o +AF /o , F h +AF fl , 
Ff 2 +AFf 2 , + A^i, and a + Aa, and consider the phase 
offset A<j>. We can expand the equation to the linear 
terms of the uncertainties Afc c , AF/ , AF/ 1? AF/ 2 , A^i, 
A0, and Aa. By neglecting the higher-order nonlinear 
terms, we can obtain 



ttlSv — —w x Sx — w v Sv - 

-kiq + wo, 
£x — —j^Sx^ q = Sx + 



2^ 

■Tf, COS 

J2 a a 



(Cl) 



where 



2tt / 

= cos 2tt — 

a \ a 

+Ak c + fc c + k 

w v = (F fl + AF fl ) /v* 

' v 



Aa 



AFt. 



w = AF fo 



In 



vi 



- sin ^2tt-- J AF /2 - F/ 2 cos (^tt— J A0. 

The characteristic equation of the linear matrix of 
Eq. (jCl|) can be expressed as: 

(s + 7/0 (s 3 + ^5 2 + w x s + fc/) = 0. 
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It can be checked from Eqs. (fT9j) and ([20]) that the 
real parts of the eigenvalues of the above equation are all 
negative. It means that there exists a stationary solution 
of the linearlization equation of Eq. (jB2|) . so does the 
original equation (|B2|) . Thus, we have q = Sx + Sx — )• 

when t —> oo. Furthermore, from &r = —jf t 5x, it can 
be verified that — > 0. Thus, we have lim^oo (5x = 
limt^oo (E (x)-x*) = 0. 
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